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Abstract — The notion of Lyapunov function plays a key role 
in design and verification of dynamical systems, as well as 
hybrid and cyber-physical systems. In this paper, to analyze 
the asymptotic stability of a dynamical system, we generalize 
standard Lyapunov functions to relaxed Lyapunov functions 
(RLFs), by considering higher order Lie derivatives of cer- 
tain functions along the system's vector field. Furthermore, 
we present a complete method to automatically discovering 
polynomial RLFs for polynomial dynamical systems (PDSs). 
Our method is complete in the sense that it is able to discover 
all polynomial RLFs by enumerating all polynomial templates 
for any PDS. 

I. Introduction 

The notion of Lyapunov function plays a very important 
role in design and verification of dynamical systems, in 
particular, in performance analysis, stability analysis and 
controller synthesis of complex dynamical and controlled 
systems [1], [2], [3]. In recent years, people realized that 
the notion is quite helpful to safety verification of hybrid 
and cyber-physical systems as well [4]. 

However, the following two issues hinder the application 
of Lyapunov functions in practice. Firstly, it is actually not 
necessary for the first-order Lie derivative of a Lyapunov 
function to be strictly negative to guarantee asymptotic 
stability, which is shown by LaSalle's Invariance Principle 
[2]. Such a condition could limit to scale up the method. 
Secondly, in general there is no effective way so far to 
find Lyapunov functions, although many methods have been 
proposed by different experts using their field expertise. 

To address the above two issues, in this paper, we first 
generalize the standard concept of Lyapunov function to 
relaxed Lyapunov function (RLF) for asymptotic stability 
analysis. Compared with the conventional definition of Lya- 
punov function, the first non-zero higher order Lie derivative 
of RLF is required to be negative, rather than its first-order 
Lie derivative. Such a relaxation extends the set of admissible 
functions that can be used to prove asymptotic stability. 

Another contribution of this paper is that we present a 
complete method to automatically discovering polynomial 
RLFs for polynomial dynamical systems (PDSs). The 
basic idea of our method is to predefine a parametric 
polynomial as a template of RLF first, and then utilize the 
Lie derivatives of the template at different orders to generate 
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constraint on the parameters, and finally solve the resulting 
constraint. Our method is complete in the sense that it is 
able to generate all polynomial RLFs by enumerating all 
polynomial templates for any PDS. 

Related Work. In [5], the same terminology "relaxed Lya- 
punov function" is used, with a different definition. 

The idea of applying higher order Lie derivatives to 
analyze asymptotic stability is not new. For example, in [6], 
[7] the authors resorted to certain linear combinations of 
higher order Lie derivatives with non-negative coefficients 
such that the combination is always negative. This method 
could be included in the framework of vector Lyapunov 
functions method [8], [9]. Our method is essentially different 
from theirs because an RLF only requires its first non-zero 
higher order Lie derivative to be negative. 

In the literature, there is a lot of work on constructing 
Lyapunov functions. For instance, in [10], [11], [12], meth- 
ods for constructing common quadratic Lyapunov functions 
for linear systems were proposed, which were generalized 
in [13] and [14] for nonlinear systems wherein the gener- 
ated Lyapunov functions are not necessarily quadratic. An- 
other useful technique is the linear matrix inequality (LMI) 
method introduced in [15] and [16], which enables us to 
utilize the results of numerical optimization for discovering 
piecewise quadratic Lyapunov functions. Based on sums-of- 
squares (SOS) decomposition and semi-definite programmng 
(SDP) [17], a method for constructing piecewise high-degree 
polynomial and piecewise non-polynomial Lyapunov func- 
tions was proposed in [18] and [19]. The SOS and SDP 
based method was also used in [20] to search for control 
Lyapunov functions for polynomial systems. In [21], the 
authors proposed a new method for computing Lypunov 
functions for polynomial systems by solving semi-algebraic 
constraint using their tool DISCOVERER [22]. Approaches 
to constructing Lyapunov functions beyond polynomials us- 
ing radial basis functions were proposed in [23], [24]. 

Our method has the following features compared to the 
related work. Firstly, it generates relaxed Lyapunov functions 
rather than conventional Lyapunov functions. Secondly, it 
is able to discover all polynomial RLFs by enumerating all 
polynomial templates for any PDS, whereas the Krasovskii's 
method [25] and Zubov's method [26] can only produce 
Lyapunov functions of special forms. Thirdly, the LMI 
method and SOS method are numerical, while our method 



is symbolic, which means it could provide a mathematically 
rigorous framework for the stability analysis of polynomial 
dynamical systems. 

Structure: The rest of this paper is organized as follows. In 
Section [III the theoretical foundations are presented. Section 
Uni shows a new criterion for asymptotic stability using the 
notion of relaxed Lyapunov functions. In Section |IV] we 
present a sound and complete method and a corresponding 
algorithm for automatically discovering polynomial RLFs on 
polynomial dynamical systems. The method is illustrated by 
an example in Section |V] Finally, we conclude this paper 
and discuss possible future work in Section [Vl] 

II. Theoretical Foundations 

In this section, we present the fundamental materials based 
on which we develop our method. 

A. Polynomial Ideal Theory 

Let K be an algebraic field, and K.[xi,X2, ■ ■ ■ ,Xn] denote 
the polynomial ring over K. Customarily, let x denote the 
?i-tuple {xi, ■ ■ ■ , Xn)- Then K[xi, X2, ■ ■ ■ , x„] can be written 
as K[x] for short, and a polynomial in K[a;i, 2:2, . . . , Xn] can 
simply be written as p(x) or p. Particularly, IK will be taken 
as the real field R in this paper, and x takes value from the 
n-dimensional Euclidean space R". 

In our method we will use polynomials with undetermined 
coefficients, called parametric polynomials or templates. 
Such polynomials are denoted by p{u, x), where u = 
(ui, W2, . . . , ut) is a t-tuple of parameters. A parametric 
polynomial p(u, x) in R[a;i, 2:2, . . . , x„] with real param- 
eters can be seen equivalently as a regular polynomial in 
R[ui,U2, ■ ■ ■ ,ut,xi,X2, ■ ■ ■ ,Xn]- Given Uq G M*, we call 
the polynomial (x) resulted by substituting Uq for u in 
p(u,x) an instantiation of p(u,x). 

The following are some fundamental results relative to 
polynomial ideals, which can be found in [27]. 

Definition 1: A subset / C ]K[x] is called an ideal iff 

(a) Oe/; 

(b) If p(x), g(x) e /, then p(x) + .g(x) G /; 

(c) If p(x) G /, then p(x)/i(x) G / for any /7,(x) G IK[x]. 
It is easy to check that if pi, • ■ • ,p,„ G K[x], then 

m 

(pi, • • • ,pm) = {^Pihi I Vi G [1, m].hi G ]K[x]} 

1=1 

is an ideal. In general, we say an ideal / is generated by 
polynomials gi, 32, ■ • • , 5fc e K[x] if J = (gi, 32, ■ • ■ , 
where all gi for i G [1, fc] are called generators of /. In fact, 
we have 

Theorem 2 (Hilbert Basis Theorem): Every ideal / C 
]K[x] has a finite generating set. That is, / = (51, 52, • • • , fffc) 
for some 51 , .92 , ■ • • , Qk G K[x] . 

From this result, it is easy to see that 

Theorem 3 (Ascending Chain Condition): For any as- 
cending chain 

/i C /2 c . . . c /„, c . . . 

of ideals in polynomial ring K[x], there must be an N such 
that for all m > N, Im = In- 



B. Dynamical Systems and Stability 

We summarize some fundamental theories of dynamical 
systems here. For details please refer to [1], [2], [3]. 

7 ) Dynamical Systems: We consider autonomous dynam- 
ical systems modeled by first-order ordinary differential 
equations 

i-f(x) , (1) 

where x G R" and f is a vector function from R" to E", 
which is also called a vector field in K". 

In this paper, we focus on special nonlinear dynamical 
systems whose vector fields are defined by polynomials. 

Definition 4 (Polynomial Dynamical System): Suppose 
f = (/i, /2, • • • , /„) in ([T]). Then ([T]i is called a polynomial 
dynamical system (PDS for short) if for every 1 < i < n, fi 
is a polynomial in R[x]. 

If f satisfies the local Lipschitz condition, then given xq G 
M", there exists a unique solution x(i) of (HJ defined on (a, b) 
with a < < 6 s.t. 

Vt G (a, b). = f (x(i)) and x(0) = xq . 

We call x(t) on [0, b) the trajectory of ([U starting from initial 

point Xq. 

Let cr(x) be a function from R" to R. Suppose both cr 
and f are differentiable in x at any order n G N. Then we 
can inductively define the Lie derivatives of a along f, i.e. 
L^'cr : R" ^ R for A: G N, as follows: 

. igcr(x) = cr(x), 

. L^a(x) = (|jL^^V(x),f(x)), for A: > 0, 
where (•, •) is the inner product of two vectors, i.e. 

( (ai, • • ■ ,a„), (61, . . . , bn) ) = I]r=i "-i^i- 

Polynomial functions are sufficiently smooth, so given 
a PDS P and a polynomial p, the vector field f of 
satisfies the local Lipschitz condition, and the higher order 
Lie derivatives of p along f are well defined and are all 
polynomials. For a parameterized polynomial p(u, x), we 
can define L^p(u, x) : R" -> R by seeing u as undeter- 
mined constants rather than variables. In the sequel we will 
impUcitly employ these facts. 

Example 5: Suppose f = {—x,y) and p{x,y) ^ x + y'^. 
Then L^p = x + y^, L\p = —x + 2y^, L|p = x + 4y^, 
L^p = —X + 8y'^ . 

2) Stability: The following are classic results of stability 
of dynamical systems in the sense of Lyapunov. 

Definition 6: A point Xe G R" is called an equilibrium 
or critical point of (HJ if f (xg) = 0. 

We assume Xe = w.l.o.g from now on. 

Definition 7: Suppose is an equilibrium of Then 

• is called Lyapunov stable if for any e > 0, there 
exists a (5 > such that if ||xo|| < S, then for the 
corresponding solution x(i) of ([TJ, ||x(f;)|[ < e for all 
t > 0. 

• is called asymptotically stable if it is Lyapunov stable 
and there exists a 6 > such that for any ||xo|| < S, 
the corresponding solution x(t) of ([T]i can be extended 
to infinity and limf_i.oo x(t) = 0. 



Lyapunov first provided a sufficient condition, using so- 
called Lyapunov function, for the Lyapunov stability as 
follows. 

Theorem 8 (Lyapunov Stability Theorem): Suppose is 
an equilibrium point of ([T]i. If there is an open set U C K" 
with E U and a continuous differentiable function V : 
U such that 

(a) y(o) = 0, 

(b) F(x) > for all x G U\{0} and 

(c) L}V{x.) < for all x G C/, 

then is a stable equilibrium of ([T]i. Moreover, if condition 
(c) is replaced by 

(c') i^F(x) < for all x e U\{0}, 

then is an asymptotically stable equilibrium of ([TJ. Such 
V is called a Lyapunov function. 

For asymptotic stability, we have Barbashin-Krasovskii- 
LaSalle (BKLS) Principle which relaxes condition (c') in 
Theorem [8] 

Theorem 9 (BKLS Principle): Suppose there exists V 
satisfying the conditions (a), (b) and (c) in Theorem |8] If 
the set = {x G M" | T/(x) = 0}nU does not contain 
any trajectory of the system besides the trivial trajectory 
x(f) = 0, then is asymptotically stable. 

Inspired by Theorem |9] we will define relaxed Lyapunov 
function (RLF for short) in the subsequent section, which 
guarantees the asymptotic stability of an equilibrium of a 
dynamical system. 

III. Relaxed Lyapunov Function 

Intuitively, a Lyapunov function requires that any trajec- 
tory starting from xq G ?7 cannot leave the region {x G 
M" I V{x) < F(xo)}. While, in the asymptotic stability 
case, the corresponding V forces any trajectory starting from 
xo G ?7 to transect the boundary {x G R" | V"(x) l/(xo)} 
towards the set {x G K" | ^(x) < V{xq)}. It is clear 
that LjV{x) < is only a sufficient condition to guarantee 
asymptotic stability. When a point x satisfies ^^^(x) ~ 0, 
the transection requirement may still be met if the first non- 
zero higher order Lie derivative of at x is negative. To 
formalize this idea, we give the following definition. 

Definition 10 (Pointwise Rank): Let N"*" be the set of 
positive natural numbers. Given sufficiently smooth function 
cr and vector filed f , the pointwise rank of a w.r.t. f is defined 
as the function ^^r.f ■ — ?> N U {oo} given by 



7<T,f(x) 



oo, if V/c G N+.L^cr(x) = 0, 
min{fc G N+ | Lf ct(x) ^ 0}, otherwise. 
Example 11: For f = {~x,y) and p{x,y) = a; + y^, 
by Example |5] we have 7pf(0,0) = co, 7pf(l,l) = 1, 
7p,f (2, 1) = 2. 

Definition 12 (Transverse Set): Given sufficiently 
smooth function a and vector field f , the transverse set of 
(J w.rt f is defined as 

Trans^,f = {x G R" | 7a,f (x) < oo A L^" '^""'' ai^x) < 0} . 

Intuitively, Trans^.f consists of those points at which the 
first non-zero high order Lie derivative of a along f is 



negative. Now we can relax condition [c] in Theorem |8] 
and get a stronger result for asymptotic stability. 

Theorem 13: Suppose is an equilibrium point of (HJ. If 
there is an open set U C R" with E U and a sufficiently 
smooth function V : U ^ M. s.t. 

(a) V{0) = 0, 

(b) V{x) > for all x G U\{0} and 

(c) X G Transy,f for all x G U\{0}, 

then is an asymptotically stable equilibrium point of 

Proof: First notice that condition (c) implies LjV{x) < 
for all X G C/\{0}. In order to show the asymptotic stability 
of 0, according to Theorem |9] it is sufficient to show that 
A^ = {x G R" I LfF(x) = 0}nU contains no nontrivial 
trajectory of the dynamical system. 

If not, let x(t), t > be such a trajectory contained in 
M other than x(t) = 0. Then LlV{x{t)) = for all t > 0. 
Noting that xo = x(0) G Transy^f, we can get the Taylor 
expansion of LjV{x{t)) at t = 0: 



LlV{x{t)) = LlV{xo) + LjVixo) ■ t + LlV{xo) ■ 



2! 



= L 



7v,f (xo) 



(7y,f(xo) - 1)! 

By Definition [121 there exists an e > s.t. Vi G 
(0, e). Lfp{x{t)) < 0, which contradicts the assumption. ■ 

Definition 14 (Relaxed Lyapunov Function): We refer to 
the function V in Theorem [13] as a relaxed Lyapunov 
function, denoted by RLF. 

In the next section, we will explore how to discover 
polynomial RLFs automatically for PDSs. 

IV. Automatically Discovering RLFs for PDSs 

Given a PDS, the process of automatically discovering 
polynomial RLFs is as follows: 

> a template, i.e. a parametric polynomial p(u, x), is 
predefined as a potential RLF; 

> the conditions for p(u, x) to be an RLF are translated 
into an equivalent formula $ of the decidable ^rsf-orc/er 
theory of reals [28]; 

> constraint on parameters u, or equivalently a set Su 
of all t-tuples subject to is obtained by applying 
quantifier elimination (QE for short. See [29], [30]) to 
$, and any instantiation of u by Ug G Su yields an RLF 

Puo(x). 

A. Computation of Transverse Set 

Correct translation of the three conditions in Theorem [T3] 
is crucial to our method. In particular, we have to show 
that for any polynomial p(x) and polynomial vector field 
f, the transverse set TranSp.f can be represented by first 
order polynomial formulas. To this end, we first give several 
theorems by exploring the properties of Lie derivatives and 
polynomial ideas. 

In what follows, given a parameterized polynomial 
p{u, x), all Lie derivatives L^p are seen as polynomials in 
R[u, x]. Besides, we will use the convention that Vie© ~ 
false and Aie0 Vi — where rji is logical formula. 



Theorem 15 (Fixed Point Theorem): Given p = p(u,x), 
if L\p e {L\p, ■ ■ ■ , L'f^p), then for all m > i, L}"p e 



V). 



Proof: We prove this fact by induction. Assume L^p e 
(LfP, ■ • ■ jLp^p) for k > i. Then there are gj S M[u,x] 
s.t. LfP = 9jLfP- the definition of Lie derivative 

it follows that 



d 



= (d^Es.Llpj) 



i=i 



d 



i-l 



= lli-d;,9,,mp + Y.g,L^'^'p 

3=1 i=i 



9 



i-l 



= '^i-g^9j,f)LiP + 9j~iLiP + 9r-lL}P- 

3 = 1 ^ i=2 

By induction hypothesis, L\p G {L^p, ■ ■ ■ , L^^^^p), so 
e (-^fP, ■ • ■ V)- By induction, the fact follows 
immediately. ■ 
Theorem 16: Given pSp(u,x), the number 

7Vp,f = min{i G N+ | L\+^p G {L\p, ■■■ , L\p)} 

is well defined and computable. 

Proof: First it is easy to show that Npj has an 
equivalent expression A'^p^f = min{i G N+ | li+i = Ii}, 
where /j = {L\p, • • • , L\p) C R[u, x]. Notice that 

/i C /2 C . . . C /fe . . . 

forms an ascending chain of ideals. By Theorem |3] iVp f 
is well-defined. Computation of Np ^ is actually an ideal 
membership problem, which can be solved by computation 
of Grobner basis [27]. ■ 
Example 17: For f = {—x,y) and p{x,y) = x + y^, by 
Example |5] we have L'^p ^ {LfP) and L^p G {L^p,Lfp), 
so Np^f = 2 . 

Theorem 18 (Rank Theorem): Suppose that p = p{u, x). 
Then for all x G R" and all Uq G K*, Tpu^.f (x) < cxd implies 
7p„o,f(x) < iVp,f. 

Proof: If the conclusion is not true, then there exist 
xo G R" and uq G R* s.t. 

^p,f < 7pu(,.f(^") < °° ■ 
By Definition [TOl xq satisfies 

ifPuo = A ■ • • A if " 'puo = A if ° puo 7^ . 

Then by Theorem [16] and [15] for all m > A'p f, we have 
i™Puo(xo) = 0. In particular, puo{xo) = 0, 



which contradicts L 



7puo,f(xo) 



Puo(xo) 7^ 0. 



Now we are able to show the computability of TranSp_f . 



Theorem 19: Given a parameterized polynomial 
p = p(u,x) and polynomial vector field f, for any Uq G R* 
and any x G R", x G TranSp^^ f if and only if Uq and x 
satisfy tpp^f, where 



'Pp. 



f ^ V ^p.f ' 



and 



(2) 



ifl, = ( /\ L',p{u, x) = 0) A L^piu, x) < . (3) 

Proof: Suppose x G TranSp^^ .f. By Definition [T2| 
X satisfies 

Ljpuo =OA---Ai;''"'"'''""'pu„ =OAL;''"«'''"puo <0 . 

(4) 

By Theorem [TSi 7p„^.f(x) < iVp f. Then it is easy to check 
that (O implies ((2| when u = Uq. 

(=>) If Uq and X satisfy ifp.f, then from Definition [T2| we 
can see that x G TranSp^^^.f holds trivially. ■ 

B. A Sound and Complete Method for Generating RLFs 



Based on the results established in Section IIV-AI we 
can give a sound and complete method for automatically 
generating polynomial RLFs on PDSs. 

Given x = {xi,X2, . . . ,a;„) G R", let ||x|| = VEI'=i ^1 
denote the Euclidean norm of x. Let i3(x, d) = {y G R" | 
||y — x|| < d} for any d > 0. Then our main result can be 
stated as follows. 

Theorem 20 (Main Result): Given a PDS x = f(x) with 
f (0) = and a parametric polynomial p^p{u, x). Let tq G 
R and Uq = U2(,, . . . , wt,,) G R*. Then p^g is an RLF 
in ^(0,7-0) if and only if 



where 



(uio,U2o, . . . ,utn,ro) G QE((/)p,f) , 



(5) 



0^^f=p(u,O) =0 , (6) 

</.2^f^^Vx.(||xf >0 A ||xf ^p(u,x) >0), (7) 

</.3_f =Vx.(||xf > A llxf < r2 ^ (^p,f) . (8) 
Proof: First, in Theorem [T3] the existence of an open 
set U is equivalent to the existence of an open set B{0, tq). 
Then according to Theorem [19] it is easy to check that (O, 
(O and dHJ are direct translations of conditions (a), (b) and 
(c) in Theorem [13] ■ 
According to Theorem [20] we can follow the three steps at 
the beginning of Section HV] to discover polynomial RLFs on 
PDSs. This method is "complete" because we can discover 
all possible polynomial RLFs by enumerating all polynomial 
templates. 



C. Implementation 

To construct (j>pj in Theorem |20l we need to compute 
Npj in advance, which is time-consuming. What is worse, 
when iVp f is a large number the resulting 0p f can be a huge 
formula, for which QE is difficult. For analysis of asymptotic 
stability, one RLF is enough. Therefore if an RLE can be 
obtained by solving constraint involving merely lower order 
Lie derivatives, there's no need to resort to higher order ones. 
Regarding this, we give an incomplete but more efficient 
implementation of Theorem |20l by constructing f and 
searching for RLEs in a stepwise manner 

Let 

j-i 

V^fS /\L^fp(u,x) = 0, fori>l , 

6i^^f = Vx.(||xf > A ||x||2 < A Vp,f 4p(u,x) < 0) 
and 

^^f^Vx.dlxf >0 A <r^ Ai';^i^L\p{u,x)<0). 

Intuitively, for x satisfying -0^ we have to impose 
constraints 0* f or ^* ^ on the i-th higher order Lie derivative 
of p along f . Now the RLE generation algorithm (RLEG for 
short) can be formally stated as follows. 

Algorithm 1: Relaxed Lyapunov Eunction Generation 



1 Input: f G m.[xi, . . . ,x„]" with f(0) 

2 peR[ui,...,ut,xi,...,x„] 

3 Output: Res C 

4 i := 1; temp := 0; L}p := {-§^p, f); 

5 i?esO := QE(0i_f A 4>l ,)- 



0, 



6 if i?es° 

7 I return 

8 else 



then 



9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 

21 return 



repeat 

temp := Res'~^ H QE(6l^ j); 
if temp ^ then 
I return temp ; 
else 

i?es* := i?es*-i nQE(^^ f); 
if Res' = then 
return 0; 

else 

z := i + 1; 
until L\p e {L\p, Lfp, 



• L\-'p); 



Remark Eormula 0^ j- and 0^ ^ in line 5 are defined in 
^ and (|7]); QE in line 5, 10 and 14 is done in a computer 
algebra tool like Redlog [29] or QEPCAD [30]; in line 20 
the loop test can be done by calling the IdealMembership 
command in Maple™ [31]. 

The idea of Algorithm [T] is: at the i-th step, we search 
for an RLE using constraint constructed from Lie derivatives 



with order no larger than i. If this fails to produce a solution, 
then we add the (i + l)-th order Lie derivative to the 
constraint. This process continues until either we succeed 
in finding a solution, or we can conclude that there is no 
RLE with the predefined template, or we get to the A'p f-th 
iteration, which means no solution exists at all. 

Correctness of the algorithm RLEG is guaranteed by the 
following theorem. 

Theorem 21: Eor Algorithm [T] we have 

1) Termination. RLFG terminates for any valid input. 

2) Soundness. If (u,r) = (ui, U2, . . . , Wt, r) £ Res, 

then Pu(x) is an RLE in S(0, r). 

3) Weak Completeness. If Res = then there does not 

exist an RLE in the form of p(u, x). 
Proof: 

(1) The loop condition is L\p ^ {L^p, L^p, . . . , V^'^^p). 
By Theorem [16] RLEG can run at most Np^f many 
iterations. 

(2) Suppose Res'-*, Res^, . . . , Res^ is the longest sequence 
generated by RLEG when it terminates. We can induc- 
tively prove that this sequence satisfies the following 
properties. 

(PI) < fc < Np,f 
(P2) Res' = QE((/.i f A 4>l, A . 
where 



for < i < fc. 



x.(|| 



> A ||xf < r^) 



(P3) QE(0i f A 02 f A 0;_f ) = 0, for 1 < i < fc, where 



;Vx. 



IxlP >0A||x||2 < 



(P4) Res = if and only if either Res^ = or fc = 
NpX, otherwise Res = QE(0i_f A 0^ f a 0^;^^). 
Suppose (u,r) G Res, then by (PI), (P4) and Q we 
can get Res C QE(0p^f). Thus (u,r) G QE(0p_f) and 
Pu(x) is an RLE according to Theorem l20l 
(3) Suppose Res = 0, then by (P4) we have either fc = Npj 
or Res^ = 0. If fc = A^p,f(> 1), then by (P3) and © 
we get QE(0p^f) = 0; if Res'^ = 0, from (PI), (P2), 
^ as well as the validity of 

Np.t k 

V<f)^((V<f) 0<k<Np,, 



we have QE(0p f ) C Res'^ = 0. So far we have proved 
Res = implies QE(0p.f) = 0. Again by applying 
Theorem |20] we get the final conclusion. 



V. Example 

We illustrate our method for RLE generation using the 
following example. 



Example 22: Consider the nonlinear dynamical system 

,2 



-X + y 



(9) 



with a unique equilibrium point 0(0, 0). We want to establish 
the asymptotic stability of O. 

First, the linearization of (|9]l at O has the coefficient matrix 



A = 







QE Vi-Vy. 



with eigenvalues —1 and 0, so none of the principles of 
stability for linear systems apply. Besides, a homogeneous 
quadratic Lyapunov function x^ + axy + hy^ for verifying 
asymptotic stability of (|9]l does not exist in R^, because 

+ y2 > ^ + axy + hy"^ > 0) 
A {2xx + ayx + axy + 2byy < 0) 

is false. However, if we try to find an RLF in using the 
simple template p = x'^+ay'^, then Algorithm[T]returns a = 1 
at the third iteration. This means ^ has an RLF x^ + y^, 
and O is asymptotically stable. 

From this example, we can see that RLFs really extend the 
class of functions that can be used for asymptotic stability 
analysis, and our method for automatically discovering RLFs 
can save us a lot of effort in finding conventional Lyapunov 
functions in some cases. 

VI. Conclusions and Future Work 

In this paper, we first generalize the notion of Lyapunov 
functions to relaxed Lyapunov functions by considering the 
higher order Lie derivatives of a smooth function along a 
vector field. The main advantage of RLF is that it provides us 
more probability of certifying asymptotic stability. We also 
propose a method for automatically discovering polynomial 
RLFs for polynomial dynamical systems. Our method is 
complete in the sense that we can enumerate all potential 
polynomial RLFs by enumerating all polynomial templates 
for a given PDS. We believe that our methodology could 
serve as a mathematically rigorous framework for the asymp- 
totic stability analysis. 

The main disadvantage of our approach is the high 
computational complexity: the complexity of the first-order 
quantifier elimination over the closed fields of reals is doubly 
exponential [32]. Currently we are considering improving the 
efficiency QE on first order polynomial formulas in special 
forms, and it will be the main focus of our future work. 
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